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ABSTRACT 


The static aeroelastic equilibrium equations for slender, straight wings are modified to incorporate the 
effects of aerodynamically-coupled formation flight. A system of equations is developed by applying trim 
constraints and is solved for component lift distribution, trim angle-of-attack, and trim aileron deflection. 
The trim values are then used to calculate the elastic twist distribution of the wing box. This system of 
equations is applied to a formation of two gliders in trimmed flight. Structural and aerodynamic 
properties are assumed for the gliders, and solutions are calculated for flexible and rigid wings in solo and 
formation flight. It is shown that formation disturbances produce greater twist in the wingtip immersed in 
the vortex than for either the opposing wingtip or the wings of a similar airplane in solo flight. Changes in 
the lift distribution, resulting from wing twist, increase the performance benefits of formation flight. A 
flexible wing in formation flight will require greater aileron deflection to achieve roll trim than a rigid 
wing. 


NOMENCLATURE 
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linear intercept 

area, meters squared 

aspect ratio 

wingspan, meters 

chord length, meters 

chord length of the aileron, meters 

lift coefficient, non-dimensional 

pitching moment coefficient, non-dimensional 

induced drag coefficient, non-dimensional 

sectional lift-curve slope, per radian 

chord length at the wing root, meters 

chord length times lift coefficient, meters 

span-wise vector of chord length times lift coefficient, meters 


| cc'j J span-wise vector of rigid-wing lift due to roll rate, meters per radian per second 

{c m } span-wise vector of pitching moment coefficient, non-dimensional 

{c m }c span-wise vector of pitching moment coefficient due to aileron deflection, per radian 

C®® J torsional influence coefficient matrix, radians per kilogram per meter 

d distance from the wing center of gravity to the aerodynamic center, meters 

e distance from the wing elastic axis to the aerodynamic center, meters 

[£] integrated influence coefficient matrix times the elastic axis, meters per kilogram 

[ /-’] integrated influence coefficient matrix times the chord squared, meters squared 

per kilogram 

g gravitational acceleration, meters per second squared 

G modulus of rigidity, kilograms per meter squared 

[G] integrated influence coefficient matrix times the center of gravity, meters per kilogram 

J torsional constant, meters to the fourth power 

k intermediate calculation 
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Subscripts 

FF 
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Superscripts 

a 

as 

e 

ib 

ob 

r 
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mass, kilograms 

span-wise vector of wing mass, kilograms 

Fourier sine series term 
load factor, non-dimensional 

roll rate, radians per second 

roll acceleration, radians per second squared 

dynamic pressure, kilograms per meter per second squared 

viscous core radius, meters 

distance, meters 

material thickness, meters 

unperturbed dependent sensitivity parameter 

perturbed dependent sensitivity parameter 

unperturbed independent sensitivity parameter 

airspeed, meters per second 

perturbed independent sensitivity parameter 

integration weighting matrix, meters 

span-wise coordinate, meters 

span-wise location of change in taper ratio, meters 

angle-of-attack, radians 

angle-of-attack at zero lift coefficient, radians 

aerodynamic operator, radians per meter 

maximum wingbox thickness in fraction of chord length 

circulation strength, meters squared per second 

aileron deflection, radians 

change in parameter 

lateral separation between wing centerlines, meters 
span-wise vector of wing twist angle, radians 
wing taper ratio 

atmospheric density, kilograms per cubic meter 
wing station 

formation flight 

inner wing section, from the wing root to the change in taper ratio 
outer wing section, from the change in taper ratio to the wingtip 
aileron deflection 


antisymmetrical 

asymmetrical 

elastic 

inboard semi-span, nearest to the leading wing 
outboard semi-span, furthest from the leading wing 
rigid 

symmetrical 
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INTRODUCTION 


For two aircraft flying in formation, certain formation geometries will produce aerodynamic 
disturbances on one or both airplanes. Generally speaking, if the airplanes are co-planar in a wingtips- 
aligned echelon formation with 3-5 wingspans of longitudinal separation, the trailing airplane will 
experience a combination of lift force and roll moment due to the wingtip vortices of the leading plane. 
By re-trimming to account for these disturbances, the trailing airplane is typically able to achieve a 
reduction in drag of 10-20% (refs. 1-3). 

Static aeroelastic deflections of the trailing wing in a formation have the potential to alter its lift 
distribution. For flexible wings, such aeroelastic effects may need to be considered when analyzing trim 
requirements or when making performance predictions. 

The majority of published literature on aerodynamically-coupled formation flight deals with 
aerodynamic modeling, automatic control system design, and to a lesser extent, performance 
measurements from wind-tunnel and flight test experiments. Very little published work is available which 
addresses methods for calculating static aeroelastic effects for a wing in formation flight. Iglesias (ref. 4) 
gives a method for optimizing the span load of wings in formation flight, but does not include elastic 
deformations. 

This paper presents a method for determining the static aeroelastic effects of formation flight 
disturbances on a slender, straight wing. Straight wings eliminate the need to consider the effects of wing 
bending on local angle-of-attack, leaving only the twisting of the wing to be accounted for in the 
aeroelastic theory. Furthermore, the effects of sweep need not be included in the aerodynamic theoiy. 
Sailplanes commonly employ slender, unswept wings and have been reported to use formation flight for 
drag reduction (refs. 5-7). Results are presented for an example sailplane. 

FORMATION FLIGHT AERODYNAMICS 

The relative geometry of two wings in an aerodynamically-coupled formation is shown in figure 1 . 
The wings are assumed to be co-planar and of sufficient longitudinal separation that there are no 
formation-induced disturbances acting on the leading wing. A coordinate system is specified such that 
span-wise coordinates are referenced to the centerline of the trailing wing. The leading wing is modeled 
as a horseshoe vortex, which consists of a bound vortex along the wing, and semi-infinite trailing vortices 
streaming from the wingtips. The trailing wing is modeled using lifting-line theory. 

The circulation, or strength per unit length, of the semi-infinite trailing vortices is assumed to be 
equal to that of the bound vortex along the wing. For longitudinal separations of three to five wingspans, 
viscous dissipation of the vortex strength can be neglected. The circulation is approximated with 
equation (1) by applying the Kutta-Joukowski theorem (ref. 8) for a wing in 1-g wings-level trimmed 
flight. 
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Figure 1. Formation geometry. 

Following standard practice for horseshoe vortices, the circulation is held constant along the span and 
set equal to the circulation at the symmetry plane of a wing with an elliptical lift distribution. This 

assumption requires that the semi-infinite trailing vortices be separated by the effective span f b . While 

this model of the wing is not exact, it does capture the fact that maximum induced drag benefit is 
achieved with some wing overlap in the span-wise direction. This real life phenomenon is due to the 
rolling up of the wing tip vortices, which is not modeled in this paper. 

For a known span-wise separation distance Ay between the wings, it is necessary to calculate the 
induced angle-of-attack at any point y along the trailing wing due to both trailing vortices of the lead 
wing. Applying the Burnham-Hallock tangential velocity profile for a vortex (ref. 9) and small-angle 
approximation, the induced angle-of-attack is shown in equation (2): 


a (y)pF ~ 


r 0 

2 n V 


Ay-fb + y 


Ay + fb + y 


(Ay-fb + y) +r} (Ay + fb + y) +r c 


( 2 ) 


The use of a viscous vortex core eliminates singularities at the vortex center. The radius r c of the 
viscous core can be chosen as 3% of the wingspan as recommended by Blake (ref. 10). 

The vortex-induced angle-of-attack distributions across a trailing wing are shown in figure 2 for three 
different formation geometries. The maximum induced upwash across the trailing wing is achieved when 
its wingtip is exactly aligned with the nearest vortex core, corresponding to the point of minimum induced 
drag. Increasing the lateral separation between the wings causes the trailing wing to miss the region of 
largest upwash, while reducing separation causes the trailing wing to encounter a significant region of 
downwash. Note that, for all three lateral separations, the angle-of-attack distributions are asymmetric 
with respect to the centerline of the trailing wing. 
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Figure 2. Vortex-induced angle-of-attack distribution. 


STATIC AEROELASTIC EQUILIBRIUM EQUATIONS 

The equations of static aeroelastic equilibrium are presented in the general matrix form used by 
Bisplinghoff (ref. 11), and are modified to include formation flight disturbances as well as to account for 
a wing in trimmed flight. The modified equilibrium equations are combined with additional trim 
constraint equations; and the resulting system is solved simultaneously for the wing lift distribution, trim 
rigid-wing angle-of-attack, and trim aileron deflection. 

General Equations Including Formation Effects 

Matrix forms of the symmetrical and antisymmetrical static aeroelastic equilibrium equations are 
given in equations (3a) and (3b) (ref. 11). These equations are augmented by adding the asymmetrical 
equation (3c) to account for formation flight effects: 
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The left-hand side of equation (3) represents the symmetrical, antisymmetrical and asymmetrical 
angle-of-attack distributions across the trailing wing. The total angle-of-attack is composed of a rigid 
wing angle (the angle-of-attack that the free-stream would make with the un-deformed wing) and a 
deformed wing twist angle. The symmetrical equation includes the effects of the rigid-wing lift, the 
associated pitching moment about the aerodynamic center, and the weight of the wing itself. 
Antisymmetrical effects include lift force and pitching moment due to aileron deflection, lift due to roll, 
and an inertial force due to roll acceleration. The asymmetrical effects are simply due to the lift obtained 
from the upwash of the leading wing. 


The aerodynamic operators 


and 


can be computed with PrandtTs lifting line theory using 


equation (4), which accounts for symmetrical and antisymmetrical effects. The asymmetrical operator 
A as must be calculated using a different theory, such as aerodynamic strip theory shown in equation 
(5), where the lift coefficient at each wing station is dependent only upon its local angle-of-attack. 
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Equilibrium Equations for Trimmed Formation Flight 

If both wings are in 1-g wings-level trimmed flight; the load factor, roll rate, and roll acceleration can 
be assumed as shown in equations (6a) and (6b): 


n 7 = 1 


(6a) 
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P = P = 0 


(6b) 


Also for trimmed flight, the symmetric aerodynamic operator can be used to express the rigid wing 
lift distribution in terms of the scalar rigid-wing angle-of-attack as shown in equation (7): 


>]{«[} >={iK <7) 


The lift due to the formation flight term can similarly be expressed using the asymmetric aerodynamic 
operator as shown in equation (8): 


={«W TO 


The resulting trimmed equations of equilibrium are shown in equations (9a), (9b), and (9c) 
respectively: 


WF =q[E]{cci} s +{!}(/ +q[F]{c m }-[G]{m w g] (9a) 
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The trailing wing-half nearest to the leading wing's vortices will be referred to as the inboard (ib) 
wing, and the wing half farthest away will be referred to as the outboard (ob) wing. The six equations that 
result when equations (9a), (9b), and (9c) are applied separately to the inboard and outboard halves of the 
trailing wing contain 10 u nk nowns. Six of these are vector unknowns consisting of the symmetrical 

{cq} s , antisymmetrical {cry ] a , and asymmetrical {cc/} flS lift distributions for both the inboard and 

outboard wing segments. The four remaining unknowns are the scalar values of trim angle-of-attack a' 
and aileron deflection 8 a for each wing half. 
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Equations (10a) and (10b) are derived from the symmetric and antisymmetric relationships between 
the two wing halves and reduce the problem to four equations and eight unknowns: 


{cc/}' ob = {cc/}' ,6 -{cc;}' (10a) 

{cc,r o1 ’ (low 


Two additional equations, (11a) and (lib), can be written to reduce the number of scalar unknowns 
from four to two. Note that the relationship in equation (lib) assumes that the deflection of the inboard 
and outboard ailerons is equal and opposite in direction. 


a r-i b =a r~ob =a r 


(Ha) 


sib sob % 
Ofl =~O a =°a 


(lib) 


The resulting system contains four equations and six unknowns. Clearly, two additional equations 
must be constructed in order to find a solution. These two equations result from the 1 -g and wings-level 
trim constraints on the trailing wing. Under the assumption that the aircraft flight path angle is very near 
zero, 1 -g trim is achieved when the component of the total lift force normal to the free-stream velocity is 
equal to the weight of the airplane. The symmetric lift terms, both rigid-wing and elastic-wing, are 
aligned according to the rigid-wing angle-of-attack and are therefore normal to the free stream. However, 
the asymmetric terms are rotated slightly forward due to the vortex-induced upwash and must be resolved 
into components relative to the free stream. Note that, a component of the associated asymmetric induced 
drag will contribute to the total aircraft lift. However, this contribution is assumed to be negligible and is 
not accounted for. Recognizing that the antisymmetrical lift cancels and the symmetrical lift is equal 
between the two wing halves, the resulting 1-g trim equation is shown as equation 12: 
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Wings-level trim is achieved when the rolling moment between the two wing halves is of equal 
magnitude and opposite sign. As in the 1-g trim case, the asymmetric lift terms are resolved into their 
components perpendicular to the free stream. The symmetric moment terms cancel and the antisymmetric 
terms sum together as shown in equation (13): 
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The system of six equations and six unknowns can now be solved using matrix inversion. The 
solution can be written as shown in equation < 1 4V 
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Wing Twist Calculations 

Once the solution to equation (14) is found, the wing torsional deflections can be computed. The total 
lift of the wing is divided into a rigid wing term and an elastic wing term as shown in equation (15): 


{cc/} = {ccf} + {ccf} (15) 


The twist angle is related to the elastic wing lift term through the appropriate aerodynamic operator as 
shown in equation (16): 


{0} = [A]{ccf} (16) 


By substituting the appropriate forms of equations (15) and (16) into the symmetrical, 
antisymmetrical, and asymmetrical trimmed equilibrium equations (9a), (9b), and (9c); equations (17a), 
(17b), and (17c) are developed to calculate the wing twist equations for a given wing segment: 

{0} s =[>]([>]-*[£]) l ^E][A s j\l}a r + q[F]{c m }-[G]{^m w g]j (17a) 


{6} a =\A a ]l\A a ]-q[E]\ l [q[E]\cc r X 8 a +q[F]{c m } § S a ) (17b) 

LJ\LJ / L 1 On a 


M"® = [ ^ ]([ A asl \-q [£]) _1 (q[E]\ A as T* {a] FF 1 (17c) 

LJ\LJ ' V, LJ 

See Appendix A for a detailed derivation of equations (17a), (17b), and (17c). 

Rigid-Wing and Solo-Flight Solutions 

For the purposes of comparison, it is desirable to obtain solutions for the same wing in solo flight as 
well as a similar but rigid wing in formation flight. For a flexible wing in solo flight, the vortex-induced 
angle-of-attack is set to zero in equation (2). Equations (14) and (17) are then solved as before. The 
solution for a rigid wing in formation flight is determined by setting the influence coefficient matrix in 
equations (3a), (3b), and (3c) to all zeros, effectively zeroing out the matrices [£], [F], and [G] . 
Equation (14) is then solved as before. 
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Performance Calculations 


An assessment of the impact of wing twist on the performance benefits of formation flight requires an 
examination of the induced drag of the trailing wing. Additional flexible wing effects on aircraft 
performance such as profile drag due to wing twist and increased aileron deflection are not addressed 
here, but may be significant. 

The induced drag for a wing is typically a function of the square of the lift coefficient. For a wing in 
formation flight, there is an additional term due to the rotation of the lift vector resulting from the vortex- 
induced angle-of-attack. The induced drag coefficient at a given wing section is calculated as shown in 
equation 18: 


C d: = ■ 


( cf + + cf* COS OCff 

kAR 


- cf sin a pp 


08 ) 


SAILPLANE EXAMPLE 

As an example, rigid wing and flexible wing trim and performance calculations will be carried out for 
a formation of two sailplanes, which are similar to the Nimbus II (Schempp-Hirth Flugzeugbau GmbH, 
Kirchheim unter Teck, Germany) open-class fiberglass sailplane designed in 1971 by Klaus Holighaus 
(ref. 12). Airfoil, performance, and geometric data for the example sailplane are listed in table 1. The 
formation is analyzed at the airspeed for best lift-to-drag ratio given in table 1 and at an altitude of 5,000 
ft. The following plane’s left inboard wingtip is assumed to be aligned with the leading plane’s right 
wingtip vortex. This is the location of maximum drag reduction, as described earlier in the Formation 
Flight Aerodynamics section. 

Other than for mass considerations when calculating trim lift, each airplane is treated as a single wing 
with no fuselage or empennage. Each wing half of the trailing plane is divided into 300 stations, and an 
integration- weighting matrix is calculated according to the Multhopp quadrature formula (ref. 11). All 
integrations are performed from the wingtip toward the wing root. 
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Table 1. Example sailplane characteristics (ref. 12). 


Characteristic 

Value 

General 


Mass: 

500 kg 

Best glide airspeed: 

105 km per hr 

Max lift-to-drag ratio: 

47.5 

Wing 


Airfoil (FX 67-K-170): 


Maximum thickness 

0.17c 

Span: 

20.3 m 

Area: 

14.41 m 2 

Chord: 


Root 

0.96 m 

Taper ratio, inner 

0.75 

Taper ratio, outer 

0.36 

Taper change point 

0.57 b 

Twist: 

0° 

Wing mass: 

230 kg 

Aileron 


Span: 


Inboard edge 

0.57 b 

Outboard edge 

0.95 b 

Chord: 

0.12m 


Structural Model 

A simple model of the sailplane wing was constructed by assuming certain structural properties. 
Bending of the wing was neglected since only the torsional influence function appears in the aeroelastic 
equilibrium equations. A rectangular single-cell wing box was assumed as shown in figure 3. 
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Figure 3. Single-cell wing box. 
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The torsional influence function can then be computed for each wing using equation (19): 
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A detailed derivation of equation (19) is contained in Appendix B. 

Constant values for the modulus of rigidity and structural wall thickness were obtained for eight plies 

of E-Glass 7781 fiberglass (ref. 13) as G = 758xlCh pounds per square inch and t = 0.0085 inches per 
ply, respectively. Values for the location of the elastic axis relative to the quarter-chord line ( e ) and to 
the chord-wise center of gravity (d ) are given in table 2 along with the wing box height (y ). As a first- 
order approximation, the distribution of wing mass along the semi-span was assumed to be proportional 
to the volume of the wing box. 


Table 2. Wing box parameters. 


Parameter 

Value 

e 

0.1250 c 

d 

-0.1250 c 

7 

0.1275 c 
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Aerodynamic Model 

Values for the necessary wing aerodynamic parameters and coefficients are summarized in table 3. 
The angle-of-attack for zero lift and pitching moment coefficient for the FX 67-K-170 airfoil were taken 
from data published by Flansman and Craig (ref. 14). The sectional lift-curve slope was approximated as 
shown in equation (20): 


„ AR 

ci =2 n 

'« AR + 2 


( 20 ) 


Table 3. Wing aerodynamic properties. 


Property 

Value 

q a 

5.8725 per radian 

a 0 

-6.33 deg 

c m 

-0.102 

Cl 8 a 

1 .25 to 1 .69 per radian 

Cm 8 a 

-0.47 to -0.49 per radian 


As recommended by Fung (ref. 15), the values for aileron lift and pitching moment coefficients were 
calculated as 80% of the theoretical values found using Glauret's approximations in equations (21) and 
(22), which are based upon the ratio of the aileron chord to the wing chord. The range of values in table 3 
for these two parameters reflects the fact that the wing chord changes along the span. 


q 8 n 



( 21 ) 


C,t % 



( 22 ) 


Results 

The symmetrical, antisymmetrical, asymmetrical, and the combined total lift distributions across the 
trailing wing are shown in figure 4. The resulting lift distribution is significantly different from an 
elliptical distribution due to the combined effects of the nearby vortex and the aileron trim. 
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Figure 4. Lift distribution for a wing in formation flight. 

Figure 5 shows a comparison of lift distributions between a flexible wing and a perfectly rigid wing, 
both in formation flight. The twisting of the flexible wing tends to reduce the lift at the root and increase 
it at the outer sections of the wing, where the ailerons are located, and for the inboard wing, where the 
vortex effects are strongest. 



Figure 5. Comparison of lift distribution for flexible and rigid wings in formation flight. 
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The component and combined total wing twist distributions for the trailing wing are shown in 
figure 6. The twist of the inboard wingtip is significantly greater than that of the outboard tip. The lift 
force generated by the vortex upwash tends to twist the inboard wingtip in a positive sense, or leading 
edge up. The trailing edge up deflection of the aileron on the inboard wing also produces a positive twist. 
Conversely, vortex upwash effects on the outboard wing are negligible, and the trailing edge down aileron 
deflection produces a negative twist on that wing half. 



Figure 6. Wing twist distribution for a wing in formation flight. 


Figure 7 shows a comparison of flexible wing twist between a wing flying in formation and a wing in 
solo flight. The wing in formation flight clearly exhibits greater twist on its inboard wing and reduced 
twist on its outboard wing. The solo wing experiences only symmetrical lift forces; therefore its twist 
distribution is symmetrical. 

A comparison of induced drag for rigid and flexible wings in both solo and formation flight is shown 
in table 4. In solo flight, the induced drag is slightly higher for the flexible wing than for the rigid wing. 
This small drag penalty is a result of the energy required to maintain the twist of the wing. 

Table 4. Induced drag. 



Rigid 

Flexible 

Solo flight 

59.6 N 

59.8 N 

Formation flight 

54.6 N 

51.0 N 
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Figure 7. Comparison of wing twist for wings in formation flight and solo flight. 

Formation flight yields an 8% reduction in induced drag for the rigid wing as compared to a 15% 
reduction for the flexible wing. This large discrepancy is attributable to the superposition of the inboard 
wing twist with the vortex-induced upwash distribution. The vortex-induced upwash causes the local lift 
vector to rotate foiward creating a component that is parallel to the free-stream and acts in a negative drag 
sense. This effect is the source of drag reduction during formation flight. In the case of a flexible wing, 
the twisting of the wing redistributes the lift away from the wing root toward the wing tips. Therefore, the 
forward-rotated lift vector near the tip is of a larger magnitude for the flexible wing than for the rigid 
wing as is the associated negative drag component. Interestingly, the deflection of the inboard wing’s 
aileron for roll trim has the opposite effect, reducing the lift near the tip along with the associated drag 
benefits. An alternate trim solution might alleviate this. 

The trim aileron deflections for flexible and rigid wings in formation flight are shown in figure 8. The 
flexible wing requires significantly greater aileron deflection to trim out the vortex-induced rolling 
moment. Reduced torsional rigidity tends to reduce aileron effectiveness. Deflection of the aileron to 
increase lift produces a downward twist of the wing leading edge, reducing the lift generated by the wing. 
The increased aileron trim deflection for a flexible wing has the undesirable effect of increasing profile 
drag on the wing, which reduces the performance benefits gained through formation flight. As an 
example, experimental results for a high lift-to-drag ratio airfoil (ref. 16) indicate a profile drag 
coefficient increase of about 0.00006 per degree of flap deflection. Table 5 shows aircraft total drag for 
rigid and flexible wings in solo flight compared to total drag in formation flight with and without aileron 
trim drag effects. Aileron trim increases the total drag by about 0.2% in the case of the rigid wing, and by 
about 0.5% in the case of the flexible wing. 
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Figure 8. Comparison of aileron trim for flexible and rigid wings in formation flight. 

Table 5. Total drag. 



Rigid 

Flexible 

Solo flight 

103.2 

103.4 

Formation flight 

98.2 

94.6 

Formation flight with aileron trim drag 

98.4 

95.1 


The scaled sensitivities of induced drag and trim aileron deflection to various aircraft parameters can 
be calculated using equation (23). Perturbed values of flexible wing induced drag coefficient ( e f / . ) and 

aileron trim deflection (S a ) were calculated for 1% changes in the wing structural and aerodynamic 
parameters listed in table 6. 


An 


<Av 


u — u V 


v — v u 


(23) 
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Table 6. Sensitivity to estimated parameters. 


V = 

u = 

c dj S a 

G*t 

0.13 

-1.9 

e 

-0.13 

0.27 

d 

1.9e-4 

-3.2e-4 

c m 

-0.03 

-3.2e-4 

a 0 

-0.30 

1.9 

% 

-0.03 

-2.7 

Cm 8 a 

0.03 

1.8 


The coefficient of induced drag is most strongly influenced by changes in the lift curve slope oq , the 
wing torsional rigidity G*t , and the location of the wing elastic axis e . For example, a 1% increase in 
ciq will reduce c ( j. by about 0.3%. A larger lift-curve slope results in a stronger lift vector rotated 

through the vortex-induced angle-of-attack. Increased torsional rigidity, and reduced spacing between the 
elastic axis and the quarter-chord line have the effect of resisting wing twist, which weakens the rotated 
lift vector and the associated reduction in induced drag. 

As to be expected, roll-axis trim requirements are strongly affected by the aileron lift and moment 
coefficients . The required aileron trim deflection is also strongly influenced by the wing torsional rigidity 
and the lift-curve slope. The wing chord- wise center of gravity and pitching moment coefficient have 
little effect on either performance benefits or lateral trim requirements. 

CONCLUSIONS 

The static aeroelastic equilibrium equations for slender, straight wings were modified to incorporate 
the effects of aerodynamically-coupled formation flight. A system of equations was developed by 
applying trim constraints and was solved for component lift distribution, trim angle-of-attack, and trim 
aileron deflection. The trim values were then used to calculate the twist distribution of the wing. The 
system of equations was applied to a formation of two gliders in trimmed flight. Structural and 
aerodynamic properties were assumed for the gliders, and the solution was calculated for flexible and 
rigid wings in solo and formation flight. 

In general, the inboard tip of a flexible wing in formation flight will tend to deflect a greater amount 
than the opposing wingtip. The inboard tip will also deflect a greater amount than the symmetrical 
deflection of a flexible wing in solo flight, while the outboard wingtip will deflect less. The outboard shift 
in lift distribution for a flexible wing has the effect of enhancing the performance benefits of formation 
flight. However, a flexible wing in formation flight will also tend to require greater roll trim input than a 
rigid wing, potentially resulting in slightly increased profile drag. 
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APPENDIX A 

DERIVATION OF EQUATIONS (17a), (17b) AND (17c) 

Begin with equations (Ala) and (Alb), based upon equations (8-83) and (8-116), respectively from 
Bisplinghoff (ref. 11). Equation (Ale) follows the general form of equation (Ala), with the only 
asymmetric effect being the vortex-induced lift. 


A {cciY =<l[ E ]{cci} S + A>< {«fj +q[F]{c m }-n z [G]{m w g} (Ala) 


A a {cciY = q[E]{cciY + A “ RL S a+\ cc l} P +l l[ F ]{ c m}§ $a-[ G \[y\{ m w} P (Alb) 
LJ L J V. L J Aj L J P J a 


A as ]{cc / } as = 9 [£]{cc / r+[A as ]{cc[}“ (Ale) 


The simplifications listed in equations (6a) and (6b) are applied along with the definitions of the 
aerodynamic operators in equations (7) and (8). The three equations can then be written as equations 
(A2a), (A2b), and (A2c), which correspond to equations (9a), (9b), and (9c): 


A S {cc/}* =q[E]{cciY + {l}« r + < l[ F ]{ c m}-[ G ]{ m wS} (A2a) 

r i ( r -i r 3 

A ° { cc lY = l[ E ]{ cc lY + A “ \ cc l\, + < i[ F ]{ c m}s s a (A2b) 

L J V L JL J °a a J 

A as ]{cciY S = 4 E ]{cciY S +{«W (A 2 c) 

Bring some terms from the left-hand side to the right-hand side: 

A ,s {cciY -q[E]{cqY -{l}«' = q[ F ]{c m }-[G]{m w g} (A3a) 
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(A3b) 


A 


{cctf-q^cctf 





+ q[ F ]{ c m} §a 




{cci} as 


q[E]{cciT S 


{ a }pF 


Factor out the {cc/} terms: 



-«[£])MMiK 


( l[ F ]{ c m}-[ G \{ m wS] 


A 


[£]){«/}“ =f[^]{»;}“ 


([ a “]-,[£]){ cc , r={a} ff 


Now expand the {cc/} terms using equation (15): 




-,[£]){«[ }’ -{IW 


^[ F ]{ c m}-[ G \{ m wS} 


([^« ] - ^ [^]) {ccf }" + ] - ^ [£]) {ecj- }" = f ] {ccf +q [F}{c m } t 




?[^]){ cc / } _ { a l FF 


(A3c) 


(A4a) 


(A4b) 


(A4c) 


(A5a) 


(A5b) 


(A5c) 
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Isolate the elastic lift coefficient distribution vectors on the left-hand side: 

{«*} =( -?[ £ ]) (- A' v + <?IXl){ccf} +{!}« r + q[ F ]{c m }-[ G ]{ m ws\ (A6a) 

{ccf} =( -?[£]) (- A a + <y[E]){ccf} + {ccf | +q[F]{c m }§ S a (A6b) 

L J \L J / L J /L J V L J L J O fl aj ^ 

{<*?}" =([A<«]- ? [£])" l f^[A“] + ,[£]){cc,'}" + {«}„,) <A6c) 

Now apply the following simplifying relations: 

V]{cc[} 5 = {1 K and {ccff = [a*] _ 1 {l}a r (A7a) 


A°]{cc[} =[a"]{cc / , '} (5 8 a and {ccf } = {ccf}^ S a 


A as {ccf} ={a} FF and {ccf} = {a} FF (A7c) 

Equations (A6a), (A6b), and (A6c) become the following: 

{ccf} A s -<?[£]} q[E] A s {1 }a r +q[F]{c m }-\G]{m w g} (A8a) 

{ccf} A a -q[E]j y[f;]{ccf}^ 8 a + q[F]{c m } § ^ S a (A8b) 

{ccf}™ =([ AflS ]-?[ £ ]) (A8c) 
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(17a) 


Finally, apply the relationship [ cc f J = [^] 1 {#} to get equations (17a), (17b), and (17c). 

{e} s = [ ]([ A 5 ] - q [£]) 1 (q[E] {1} a' +q[F] {c m } - [G] {± m w g \ ) 

{0} a =[^ a l([^ a l-?[£]r 1 U[^]W}“ S a +q[F]{c m } s S a ) (17b) 

L J\L J ' V, L a J 

{0r=[A‘“]([A‘“]- 9 [£])" 1 ^[£][A as ]" 1 {a} FF j (17c) 
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APPENDIX B 

DERIVATION OF EQUATION 19 

The torsional influence function can be computed for each wing from equation 2-120 of Bisplinghoff 
(ref. 11): 




bl 2 


dy 

GJ 


(Bl) 


In order to evaluate the torsional influence function integral equation (Bl) directly, the torsional 
rigidity GJ must be written in terms of the span-wise coordinate y . For a single cell box, the torsional 
rigidity is computed from equation 2-125 of Bisplinghoff (ref. 11): 


GJ = G 


4 A- 


ds 

t 


(B2) 


From the wing box geometry shown in figure 4, the dimensions of the wing box are: 


length = \ c 
height = y c 


(B3) 


Under the assumption that the material thickness is constant, the integral in the denominator of 
equation (B2) for a rectangular wing box is simply the perimeter divided by the material thickness as 
shown in equation (B4): 


d) 


ds 

t 



(B4) 


The area of the wing box is easily found as: 


T = 



(B5) 


The rigidity equation can now be written in terms of chord length by substituting equations (B4) and 
(B5) into equation (B2): 
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GJ = 


4 tG 

{yc)\ 

(l: 

0 2 

2 

(yc+ 




-> GJ = ^ 
2y + l 


(B6) 


Chord length is related to the span-wise coordinate y through the taper ratio. The example sailplane 
has a two-piece wing with different taper ratios for the inner and outer sections. Let y(AA) indicate the 
point along the wing where the taper ratio changes and c r be the chord length at the wing root. The inner 
and outer taper ratios are then defined as follows: 


k = 


c(y) 


,a,«JA 0 .dz) 


'y=AX 


y=\b 


(B7) 


The chord length of each wing section varies linearly with the span, according to equation (B8): 


f dc ^ 


c(v),- = (^— j v + a; , y<y(Av) 
c (y)o = [ C G~\ y +a o > y>y( A y) 


\ d yj 


(B8) 


Evaluating the slope and intercept for the inner section: 


r dc} c Wly=o c (- v )ly=v(AT) 


\dyj 


0 -y(Ah) 


Cf C)~ Xj 

y(AT) 


(B9) 


a i = c (y)Li=c 


— > a i = c r 


(BIO) 


The chord length for the inner wing half, then, is: 


c (>’),• = c r 


1- 


v y(AA) y 


0 < y < y(AX) 


(Bll) 
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The outer wing half equation is similarly developed by evaluating equation (B-8) at the y{AX) station 
and at the wing tip, y = \b : 


c (y ) 0 

c (y ) 0 


'y=y{Ay) 


' dc 
ydy 


\ y=\b 


( dc A 


dy 


' o 


\ 

y(AX) + a 0 - c r X t 
) o 



c r X () 


Solving these two equations simultaneously yields: 


(B 12) 


(dc} c r (Xj-X 0 ) \bXi-X 0 y(AX) 

— = f 2 - , and a n = c r — 

V dyJ 0 y(AX)-^b ^b-y(AX) 


(B 13) 


The chord length of the outer wing half, then, is: 


c (y)o= c r 


2 bX{ ~ X 0 y(AX) 

y i b -y( A V , 


f 




hj - h Q 
hb-y(AX)^ 


y 


y{AX)<y<\_b 


(B14) 


The constants are defined: 


k _ tG Y 2 Cr 

2y 2 + 1 


(B15a) 


k; = 


i -hj 
y(AX) 


k 1 = 2 

K-O — 


\ bXj - X 0 y(AX) 
\b - y{AX) 


2 = h, - X 0 

Kq — 


hb-y(AX) 


(B 1 5b) 


(B15c) 


(B15d) 
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Equations (B 11) and (B14) are then written as: 


iy)= 


c r [\-kiy\ , 0 <y<y(M) 

y{AX)<y<\_b 


K o k a y 


(B16) 


Substituting equations (B6) and (B16) into equation (Bl) and integrating piecewise over the inner and 
outer wing sections gives the torsional influence function C . The inner section integral is evaluated as: 



r d y - 1 

r d y - 1 ft 


)'°J' 2 c 3 l 
2y + 1 1 

i 

) 

i i 

Q 

+ ^ 

1 1 

~ k J v 


3 d y 


let 


u = 1 - kj y 
du = -kjdy 


then 


c ee 


-l 

kkj 


y 

J u ^ du 
0 


cf e 


-1 

kk: 


-2 


JO 


c ee 


l ~ 

2kki . 


(i -by) 


-2 



(B17) 
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The outer section integral is evaluated similarly: 


let 


C 


ee 


I 


dy 


1 


dy 


y(AX )^ c 3 y (AX) 
2y + l ° 


tGy~ 
2y + 1 


c r k 0 


o - *£y) 


-i 3 k 


\ j [ k o - k ly) dy 


_v( 4 A) 


= J 


dy 


- / 


dy 


y(AX)^ 'K_ c l y(AX) 
2y + l 


tGy" 
2y + 1 


c r I k 0 


o - k h) 


3 k 


= T J \ k 'o ~ k ly) dy 


v(zlA) 


let 


then 


U = kg — kgV 
2 

du — —kgdy 


r 90 _ 
c o 2 

kk O y(AX) 


ee _ -i 


J u du — ^ — 


kkX 


-2 


-2 


r ee _ 


2 kkX 


y 

M a ^) 

-2 


— » 


[ k o~ k ly] -[ k \,-K 


-2 


(B18) 


Equation ( 1 9), the piecewise expression for the torsional influence function, follows from equations 
(B 17) and (B18): 


c ee = i 

2 kki 


C 96 =C de {y {AX)) - 


,(!-M 


--1 


0 < y < y{AX) 


2 kkx 


[kl-k 2 0 y] (/ k l 0 -k 2 0 y{AX )) 


y{AX)<y<\b 


(B19) 


30 



REPORT DOCUMENTATION PAGE 


Form Approved 
OMB No. 0704-0188 


The public reporting burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing instructions, searching existing 
data sources, gathering and maintaining the data needed, and completing and reviewing the collection of information. Send comments regarding this burden estimate or 
any other aspect of this collection of information, including suggestions for reducing this burden, to Department of Defense, Washington Headquarters Services, 
Directorate for Information Operations and Reports (0704-0188), 1215 Jefferson Davis Highway, Suite 1204, Arlington, VA 22202-4302. Respondents should be aware 
that notwithstanding any other provision of law, no person shall be subject to any penalty for failing to comply with a collection of information if it does not display a 
currently valid OMB control number. 

PLEASE DQ NQT RETURN YQUR FORM TQ THE ABOVE ADDRESS ¥ 


1. REPORT DATE (DD-MM-YYYY) 2. REPORT TYPE 

01-08-2009 Technical Memorandum 

3. DATES COVERED (From - To) 

4. TITLE AND SUBTITLE 

Static Aeroelastic Effects of Formation Flight for Slender Unswept Wings 

5a. CONTRACT NUMBER 

5b. GRANT NUMBER 

5c. PROGRAM ELEMENT NUMBER 

6. AUTHOR(S) 

Curtis E. Hanson 

5d. PROJECT NUMBER 

5e. TASK NUMBER 

5f. WORK UNIT NUMBER 

7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 

NASA Dryden Flight Research Center 
P.O. Box 273 

Edwards, California 93523-0273 

8. PERFORMING ORGANIZATION 
REPORT NUMBER 

H-2980 

9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESS(ES) 

National Aeronautics and Space Administraction 
Washington, DC 20546-001 

10. SPONSORING/MONITOR'S ACRONYM(S) 

NASA 

11. SPONSORING/MONITORING 
REPORT NUMBER 

NASA/TM-2009-2 14649 


12. DISTRIBUTION/AVAILABILITY STATEMENT 

Unclassified — Unlimited 

Subject Category 02, 05, 39 Availability: NASA CASI (301) 621-0390 Distribution: Standard 


13. SUPPLEMENTARY NOTES 

Hanson, NASA Dryden Flight Research Center. 


14. ABSTRACT 

The static aeroelastic equilibrium equations for slender, straight wings are modified to incorporate the effects of aerodynamically- 
coupled formation flight. A system of equations is developed by applying trim constraints and is solved for component lift 
distribution, trim angle-of-attack, and trim aileron deflection. The trim values are then used to calculate the elastic twist distribution 
of the wing box. This system of equations is applied to a formation of two gliders in trimmed flight. Structural and aerodynamic 
properties are assumed for the gliders, and solutions are calculated for flexible and rigid wings in solo and formation flight. It is 
shown for a sample application of two gliders in formation flight, that formation disturbances produce greater twist in the wingtip 
immersed in the vortex than for either the opposing wingtip or the wings of a similar airplane in solo flight. Changes in the lift 
distribution, resulting from wing twist, increase the performance benefits of formation flight. A flexible wing in fonnation flight will 
require greater aileron deflection to achieve roll trim than a rigid wing. 


15. SUBJECT TERMS 

Aeroelastic wing twist, Drag reduction. Formation flight, Sailplane, Wingtip vortices 


16. SECURITY CLASSIFICATION OF: 

17. LIMITATION OF 
ABSTRACT 

18. NUMBER 
OF 

19b. NAME OF RESPONSIBLE PERSON 

a. REPORT 

b. ABSTRACT 

c. THIS PAGE 

PAGES 

STI Help Desk (443) 757-5802 





19b. TELEPHONE NUMBER (Include area code) 

u 

u 

u 

uu 

35 

(443) 757-5802 


Standard Form 298 (Rev. 8-98) 

Prescribed by ANSI Std. Z39-18 

































